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Abstract. We show the results ol two-dimensional, azimuthally symmetric simula- 
tions of advection dominated fluid flows around black holes. We use the SPH method, 
and the a-prescription for the viscosity, including all terms in the stress tensor. 
We have performed calculations with 0.001 < a < 0.2 and find strong circulation 
patterns in the flow, on spatial scales that depend on the magnitude of the viscosity. 
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1. Introduction, method and initial conditions 

The astro-physical relevance of accretion flows in which radiative cooling 
is inefficient has been recently emphasized (Narayan and Yi, 1994). 
These flows show complicated hydrodynamical behavior that is very 
different from that obtained with the thin disk approach (Shakura and 
Sunyaev, 1973) for efficient cooling. 

We use the Smooth Particle Hydrodynamics (SPH) method (Mon- 
aghan, 1992), in azimuthal symmetry. The equation of state is P = 
pu(7 — 1) with 7 = 4/3. Viscosity is computed using an a prescription 
(/i = apc^/fifc), including all the terms from the viscous stress tensor. 
The black hole is modeled as a Newtonian point mass with an absorbing 
boundary at r in = 1.3r G = 1.2,{2GM BH / c 2 ) . We use two kinds of 
initial conditions: 1) equilibrium hydrostatic tori with constant specific 
angular momentum; 2) matter injection at a steady rate at a radius 
r = 300rc. We assume no radiative energy losses. 



2. Results 

For initially hydrostatic tori, the presence of a hot, high-entropy bubble 
near the equatorial plane is easily seen (Figure 1). This moves outward 
in the disk, on top of a background accretion flow. There is a large-scale 
circulation pattern, on a scale z — r. For a lower viscosity a = 0.001 
there are small-scale vortices in the disk, and the overall behavior is 
more turbulent. A background accretion flow is also present. 
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For the case of matter injection, the flow structure is governed by the 
magnitude of the viscosity, with small-scale eddies at low values of a 
and large-scale circulations when the viscosity is high. The fluid shows 
significant oscillations in the vertical direction as it moves radially 
inward. 

With these simulations we have tested our code with previously 
published results obtained with a different numerical scheme (Igumen- 
shchev et al., 1996; Igumenshchev and Abramowicz, 1999), and we have 
found excellent agreement both in a qualitative and quantitative man- 
ner. In future work we will extend our results to study the parameter 
space (a, 7) and increase the resolution. 
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Figure 1. Momentum vector field pv for a = 0.01. 
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